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ABSTRACT 

An  important  issue  presented  here  is  the  existence  of  two  different  types  of  correlation  between  time  rate  of 
change  of  sea  surface  temperature  (SST)  dTJdt  and  the  ocean  mixed-layer  (OML)  depth  hw,  depending  on  the 
ocean  surface  conditions.  In  the  tropics,  the  net  heat  flux  (including  radiation  budget)  at  the  ocean  surface  is 
generally  downward  and  dampens  the  OML  turbulent  kinetic  energy  (TKE).  As  the  ocean  surface  is  under 
weak  wind  forcing,  the  OML  shallows  to  the  Monin-Obukhov  length  scale  due  to  insufficient  wind-generated 
TKE  to  entrain  the  deep  cool  water  into  the  OML.  Net  heat  gain  at  the  ocean  surface  and  very  little  heat  loss 
at  the  OML  base  leads  to  a  warming  OML.  The  thinner  the  OML,  the  warmer  the  OML.  Therefore,  dTJdt 
and  hw  are  negatively  correlated;  however,  as  the  ocean  surface  is  under  strong  wind  forcing,  the  OML  deepens 
due  to  sufficient  wind-generated  TKE  to  entrain  the  deep  cool  water  into  OML.  The  excessive  heat  loss  at  the 
OML  base  (entrainment  heat  flux)  over  net  heat  gain  at  the  ocean  surface  leads  to  a  cooling  OML.  The  thicker 
the  OML,  the  warmer  the  OML.  Therefore,  dTJdt  and  h„  are  positively  correlated.  The  shift  of  positive-only 
correlation  between  dTJdt  and  hw  (as  conventionally  thought  to  be),  to  two-way  correlation  leads  to  a  new 
theory  about  the  interannual  tropical  ocean-atmosphere  interaction. 

To  investigate  the  impact  of  this  new  idea  on  the  interannual  tropical  ocean-atmosphere  interaction,  an  Air- 
Ocean  Surface  Heat  Exchange  (AOSHE)  Model,  which  is  a  coupled  system  consisting  of  the  Wind-Induced 
Surface  Heat  Exchange  (WISHE)  Model  (Yano  and  Emanuel)  and  the  Ocean  Mixed-Layer  (OML)  Model 
(Chu  and  Garwood)  has  been  developed.  The  major  advantage  of  the  AOSHE  model  is  simple  but  contains 
realistic  thermodynamics  in  both  ocean  and  atmosphere.  The  SST,  predicted  by  the  OML  model,  largely  affects 
the  surface  evaporation  rate,  which  in  turn  impacts  on  the  atmosphere.  On  the  other  hand,  the  surface  wind 
stress  and  buoyancy  flux,  predicted  by  the  WISHE  model,  drives  the  ocean  mixed  layer.  Two  different  types  of 
low-frequency  (interannual)  modes  are  found  in  the  AOSHE  model  for  different  ocean  surface  conditions.  For 
weak  surface  wind  forcing,  the  low-frequency  mode  (called  AOSHE  W-mode)  propagates  westward  with  phase 
speed  around  20  cm  s-1,  and  maximum  growth  rate  2.2  yr_l.  The  growth  rate  increases  as  the  wavenumber 
decreases.  The  most  unstable  mode  is  at  the  lowest  wavenumber  with  periods  of  2-10  years.  For  strong  surface 
wind  forcing,  however,  the  low-frequency  mode  (called  AOSHE  S-mode)  propagates  eastward  with  phase  speed 
10-100  cm  s“‘  and  period  1-3  years.  This  low-frequency  mode  has  strong  wavelength  selection.  The  maximum 
growth  rate  (nearly  0.5  yr~')  appears  at  the  3-5  zonal  wavenumber. 


1.  Introduction 

The  phenomenon  of  El  Nino  and  Southern  Oscil¬ 
lation  (ENSO)  exhibits  a  time  scale  of  2-9  years.  Two 
complementary  phases  (La  Nina  and  El  Nino)  of  the 
Southern  Oscillation  are  defined  by  Philander  ( 1985). 
Several  important  mechanisms  (virtually  wave-type) 
of  air-ocean  interaction  for  the  ENSO  phenomenon 
have  been  developed  in  the  past  decade:  the  slow  prop¬ 
agation  of  oceanic  Rossby  and  Kelvin  waves  across  the 
Pacific  Ocean  (e.g.,  McCreary  1983;  Battisti  1988;  Bat- 
tisti  and  Hirst  1989;  Cane  et  al.  1990;  Graham  and 
White  1991);  the  generation  of  two  states — El  Nino 
and  La  Nina  (e.g.,  McCreary  and  Anderson  1984); 
and  air-ocean  coupled  instability  (e.g.,  Hirst  1986; 
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Rennick  and  Haney  1986;  Yamagata  1985;  Zebiak  and 
Cane  1987;  Kuo  1989).  Despite  the  recent  debate  be¬ 
tween  Graham  and  White  ( 1991 )  and  Battisti  ( 1991 ) 
on  the  role  of  off-equatorial  oceanic  Rossby  waves  dur¬ 
ing  ENSO,  and  despite  some  speculations  on  ENSO 
events  triggered  by  high-frequency  (30-60  day)  oscil¬ 
lations  (e.g.,  Lau  and  Chan  1986;  Lau  and  Shen  1988; 
and  Hirst  and  Lau  1990),  the  ENSO  oscillations  are 
usually  thought  to  be  driven  by  delayed  negative  feed¬ 
back  due  to  reflected  oceanic  Rossby  wave  activity, 
aspects  which  can  be  traced  through  Wyrtki  (1975), 
McCreary  ( 1983),  Cane  and  Zebiak  ( 1985),  and  Bat¬ 
tisti  ( 1988,  1989)  among  others.  These  models,  how¬ 
ever,  involve  two  major  assumptions:  (a)  SST  pertur¬ 
bation  or  the  time  rate  of  change  of  SST  perturbation 
is  positively  correlated  to  the  perturbation  of  upper- 
ocean  thickness  (the  thicker  the  upper-ocean  layer,  the 
higher  the  SST);  (b)  surface  winds  over  the  eastern 
Pacific  are  switched  from  equatorial  easterlies  repre- 
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senting  Walker  circulation  when  the  eastern  Pacific  is 
cool  (La  Nina)  to  extra-equatorial  (off-equatorial) 
easterlies  representing  enhanced  Hadley  circulation 
when  the  eastern  Pacific  is  warm  (El  Nino).  This  was 
presented  by  Bjerknes  ( 1 966,  1 969 )  and  later  pursued 
in  numerical  models  by  McCreary  ( 1983),  Zebiak  and 
Cane  (1987),  and  Battisti  ( 1988,  1989).  These  are  the 
wave-type  ENSO  theories. 

Recently,  a  SST-type  ENSO  theory  was  presented 
by  Neelin  ( 1991 ).  The  new  idea  brought  to  the  com¬ 
munity  is  the  separation  of  the  slow  SST  mode  and 
the  (relatively  fast)  modified  oceanic  wave  modes. 
Within  the  SST  mode,  the  ocean  dynamics  tends  to 
come  into  balance  quickly  on  the  interannual  time 
scales  over  which  SST  is  evolving.  Neelin  (1991)  further 
argues  that  it  is  the  slow  SST  mode  rather  than  the 
time  delay  associated  with  oceanic  wave  propagation 
across  the  basin  that  is  essential  to  the  existence  of 
interannual  coupled  oscillations.  If  Neelin’s  conclusion 
is  correct,  the  thermodynamics  in  ocean  and  atmo¬ 
sphere  should  play  a  key  role  in  determining  the  SST 
mode.  In  Neelin’s  model,  the  ocean  component  is  based 
on  a  modified  shallow-water  model  with  a  fixed-depth 
mixed  layer.  The  SST  equation  contains  advection 
(both  horizontal  and  vertical)  and  a  Newtonian  cooling 
term  representing  all  one-dimensional  vertical  pro¬ 
cesses:  mixing,  sea  surface  sensible  and  latent  heat 
fluxes,  long-  and  shortwave  radiation.  The  atmospheric 
part  is  a  functional  relationship  between  the  surface 
wind  stress  and  SST.  In  the  real  ocean,  the  OML  depth 
has  large  temporal  and  spatial  variations,  and  the  one¬ 
dimensional  vertical  processes  are  too  complicated  to 
be  interpreted  as  a  Newtonian  cooling  term.  The  at¬ 
mospheric  thermodynamics  is  also  too  simple  to  be 
realistic  for  the  problems  for  which  thermodynamical 
processes  are  important.  Therefore,  the  question  may 
arise:  What  is  the  behavior  of  the  SST  mode  in  a  more 
thermodynamically  realistic  coupled  system? 

For  the  atmospheric  thermodynamics  in  most  ENSO 
models  it  is  simply  assumed  that  evaporation  from  the 
ocean  surface  is  immediately  reflected  as  heating  of  the 
atmosphere,  which  implies  upward  velocity  everywhere 
in  the  atmosphere.  In  investigating  the  linkage  between 
the  40-50-day  oscillation  and  ENSO  (Lau  and  Shen 
1988),  the  atmospheric  thermodynamics  is  improved 
by  expressing  precipitation  as  the  sum  of  surface  evap¬ 
oration  and  moisture  convergence.  Their  study  shows 
the  importance  of  moist  processes  1 )  in  generation  of 
eastward-propagating  unstable  modes  on  an  intrasea- 
sonal  time  scale  (30-50  days)  and  2)  possibly  to  onset 
and  evolution  of  ENSO-like  events  in  the  tropical  at¬ 
mosphere-ocean  system.  A  coupled  air-ocean  system 
with  an  atmospheric  part  from  the  Lau  and  Shen 
(1988)  model  and  an  oceanic  part  from  the  Hirst  ( 1 986 ) 
model  was  developed  (Hirst  and  Lau  1 990)  to  determine 
the  relationship  between  the  intraseasonal  and  inter¬ 
annual  oscillations.  The  Hirst  and  Lau  model  results 
emphasize  the  importance  of  including  prognostic  at¬ 


mospheric  equations  in  simple,  coupled  ocean-atmo¬ 
sphere  models.  The  atmospheric  part  of  the  Lau  and 
Shen  ( 1988)  and  Hirst  and  Lau  ( 1990)  models  and  any 
other  simple,  coupled  models  still  imply  an  assumption 
that  vertical  velocity  in  the  atmosphere  is  upward  ev¬ 
erywhere,  such  that  moisture  evaporated  from  the  ocean 
surface  can  directly  warm  the  atmosphere.  The  oceanic 
part  of  the  Hirst  (1986)  and  Hirst  and  Lau  ( 1 990 )  mod¬ 
els  contains  an  OML  in  one  of  the  two  regimes  (to  be 
discussed  in  detail  in  section  2). 

Schopf  and  Suarez  (1988)  proposed  a  coupled  dry 
atmosphere  and  embedded  mixed  layer-ocean  circu¬ 
lation  model  for  simulating  ENSO-like  variability.  The 
model  predicts  some  incorrect  results;  that  is,  the  ocean 
mixed  layer  is  cooler  than  surface  air  over  most  of  the 
tropical  ocean  basin  ( model  climatology ) .  Several  de¬ 
cades  ago,  meteorologists  were  already  aware  of  the 
fact  that  the  largest  part  of  the  tropical  oceans  appears 
to  be  slightly  warmer — near  0.5°C  on  the  average — 
than  the  surface  air,  both  in  the  annual  mean  and  sea¬ 
sonal  breakdowns  (Riehl  1954).  Recent  analysis 
( Hsiung  1986)  from  a  comprehensive  dataset  covering 
the  global  ocean  surface  for  a  period  from  1946  to  1979 
provided  by  Fleet  Numerical  Oceanographic  Center 
also  shows  similar  results:  the  tropical  ocean  is  warmer 
(0. 1  °-0.9° )  than  surface  air.  Two  flaws  in  the  Schopf- 
Suarez  model  (1988)  cause  this  discrepancy:  1)  the 
thermodynamics  is  improperly  parameterized  so  that 
the  atmosphere  is  thermally  driven  by  both  radiation 
(Newtonian  “ cooling ”  form)  and  air-ocean  interfacial 
heat  flux,  but  the  ocean  is  thermally  forced  only  by  the 
interfacial  heat  flux;  and  2)  wind  stress  is  used  in  two 
contradictory  ways:  surface  shear  production  S{  wrong 
dimension  used  by  Schopf  and  Suarez  1988)  is  used 
as  a  constant  in  ocean  surface  TKE  production  and  as 
a  variable  (proportional  to  the  lower-layer  wind  veloc¬ 
ity)  in  ocean  surface  current  generation.  Furthermore, 
even  the  Coriolis  force  has  the  wrong  sign  in  the  ocean 
model.  In  fact,  observations  show  that  incoming  solar 
radiation  and  latent  heat  flux  are  the  two  dominant 
components  controlling  net  surface  energy  fluxes 
(Hsiung  1986),  that  is,  the  tropical  ocean  is  mainly 
heated  by  solar  radiation  and  cooled  by  evaporation. 
Heat  gain  in  the  OML  is  not  the  same  as  the  heat  loss 
from  the  atmosphere.  Absence  of  solar  radiation  in  the 
ocean  part  of  the  Schopf-Suarez  model  leads  to  this 
unrealistic  result  that  the  OML  is  cooler  than  the  sur¬ 
face  air  over  most  of  the  tropical  ocean  basin.  More¬ 
over,  the  assumed  constant  ocean  surface  turbulent 
production  by  wind  mixing  ( 1 .25  X  10~6  m3  s-3),  cor¬ 
responding  to  a  rather  strong  surface  wind  (~10 
m  s~'),  excludes  the  effect  of  weak  surface  winds  on 
OML  dynamics,  and  the  effects  of  surface  wind  vari¬ 
ation  on  thermodynamics  of  the  atmosphere  (wind- 
evaporation  feedback,  e.g.,  Emanuel  1987)  and  oceans 
(wind-OML  feedback). 

Therefore,  to  investigate  the  SST  mode,  more  real¬ 
istic  thermodynamics  in  both  ocean  and  atmosphere 
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are  needed.  The  SST  mode  is  further  investigated  in 
this  research  by  a  coupled  model  with  more  realistic 
thermodynamics  in  both  atmosphere  and  ocean:  the 
Air-Ocean  Surface  Heat  Exchange  (AOSHE)  Model, 
which  is  a  coupled  system  consisting  of  the  Wind-In¬ 
duced  Surface  Heat  Exchange  ( WISHE)  Model  ( Yano 
and  Emanuel  1991;  hereafter  YE91),  and  the  OML 
model  (Chu  et  al.  1990;  hereafter  CGM90;  Chu  and 
Garwood  1991;  hereafter  CG91 ). 

2.  AOSHE  Model 

Figure  1  presents  a  schematic  view  of  the  coupled 
WISHE  and  OML  model  (the  AOSHE  Model).  Above 
the  air-ocean  interface,  it  is  the  ordinary  WISHE 
Model  depicted  by  E87  and  YE9 1 .  Below  the  air-ocean 
interface,  it  is  the  OML  Model  described  by  CGM90 
and  CG91,  which  is  a  modified  Kraus  and  Turner 
( 1967)  model.  The  fluxes  across  the  interface  are  no 
longer  fixed  and  drive  motions  in  both  atmosphere  and 
ocean. 

a.  An  OML  switcher 

Arguments  are  cast  in  terms  of  simple  mixed-layer 
models  (Chu  et  al.  1990;  Chu  and  Garwood  1991), 
where  it  is  assumed  that  the  temperature  and  velocity 
are  uniform  over  some  depth  hw,  the  thickness  of  the 
OML,  and  that  the  penetration  depth  of  solar  radiation 
is  much  smaller  than  hw.  With  these  assumptions,  one 
can  write 


dTs 

— -  +  uw-VTs  = 
at 


Qo  ~  Lw(J\)Q-h 

PwOCpwhw 


(1) 


where  Ts  is  the  SST,  u„,  is  the  horizontal  velocity,  pw0 
is  the  characteristic  seawater  density,  cpw  is  the  seawater 
specific  heat  under  constant  pressure,  Q()  is  the  net  sur¬ 
face  heat  flux  with  downward  being  positive,  AW(II)  is 
the  Heaviside  step  function  of  II,  and  II  is  defined  by 

II  =  C\UW*  C2otghwQo/  pwoCpw>  (2) 

denoting  surface  forcing  condition.  Here  C,  and  C2  are 
tuning  coefficients,  a  is  the  thermal  expansion  coeffi¬ 
cient, 


uw*  =  (Cz>pao/Pwo)l/2|VA|  (3) 


is  the  water  surface  friction  velocity,  and  |  V*  |  is  the 
subcloud-layer  wind  speed;  Q-h  is  the  entrainment  heat 
flux  at  the  base  of  the  OML: 


Q-h  , 

- =  we(Ts  -  T-h),  (4) 

PwO^pw 

where  we  and  T-h  are  the  entrainment  velocity  and  the 
temperature  at  the  base  of  the  OML. 

Entrainment  velocity  can  also  be  parameterized  in 
terms  of  OML  turbulent  kinetic  energy  (TKE)  balance. 
Since  salinity  effects  are  not  considered  at  present, 
sources  and  sinks  of  TKE  at  ocean  surface  are  wind 


work  (proportional  to  wind  speed  cubed)  and  buoyant 
damping  (forcing)  due  to  surface  warming  (cooling). 
This  is  given  by 


gnwOCyTs  -*—/?) 

The  step  function  Aw  acts  as  an  OML  switcher. 

The  value  of  Aw  is  1  for  II  >  0,  which  is  usually 
associated  with  strong  surface  wind  forcing  in  tropics 
(Qo  >  0).  The  heat  equation  ( 1 )  for  the  entrainment 
regime  becomes 


dTs  ,  Q-h  ~  Qo 

—  +  u„,  •  VTS  = - — 

ut  PwQCpw*lw 


(6a) 


In  this  case,  the  wind-generated  TKE  is  strong  enough 
to  entrain  the  deep  cooler  water  into  the  OML.  The 
heat  loss  at  the  OML  base,  that  is,  the  entrainment 
heat  flux  ( Q-h),  is  larger  than  the  heat  gain  at  the  ocean 
surface:  Q-h  >  Qo  (Miller  1976),  which  makes  the 
OML  a  cooling  layer.  As  shown  in  (6a),  an  increase 
in  the  OML  thickness  hw  leads  to  an  increase  of  BTJ 
dt,  which  implies  a  positive  correlation  between  dTJ 
dt  and  hw.  The  entrainment  rate  is  determined  by  (5) 
and  is  substituted  into  the  following  equation  to  prog¬ 
nosticate  the  mixed-layer  depth  hw : 


dhw 

— -  =  we-  W-h. 


(7a) 


Here  w-h  is  the  vertical  velocity  at  the  mixed-layer  base 
caused  by  Ekman  pumping  or  by  the  ocean  interior 
dynamical  processes.  Another  way  to  view  (7a)  is  that 
the  entrainment  velocity  we  is  the  summation  of  the 
vertical  velocity  at  the  OML  base  w-h  and  the  time 
rate  change  of  OML  thickness  dhw/dt.  The  vertical 
thermal  advection  in  OML  base  becomes  part  of  the 
heat  flux  Q-h  through  (4). 

The  value  of  A„,  is  0  for  II  <  0,  which  is  usually 
associated  with  weak  surface  wind  forcing  in  the  tropics. 
In  this  case,  the  wind  generated  TKE  is  so  weak  that 
it  is  unable  to  entrain  the  deep  cooler  water  into  the 
OML  ( Q-h  =  0) .  The  OML  then  shoals  to  the  Monin- 
Obukhov  length  scale.  The  heat  equation  ( 1 )  for  the 
shallowing  regime  becomes 


dT, 

dt 


+  uw-VTs 


Qo 

PwOCpwhw 


(6b) 


The  mixed-layer  depth  is  calculated  from  a  balance  of 
the  remaining  terms  in  (5),  and  it  equals  the  Monin- 
Obukhov  length  scale  /„,: 


I  _  0\  u  w* 

C2agQ0/  PwoCpw 


(7b) 


Absence  of  vertical  advection  in  heat  equation  (6b) 
for  the  weak  surface  wind-forcing  case  ( Aw  =  0)  is  due 
to  the  small  cooling  rate  caused  by  vertical  advection. 
The  Ekman  layer  thickness,  hE  =  {2K/  |/|)1/2,  in- 
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creases  equatorward  and  tends  to  infinity  at  the  equa¬ 
tor.  Here,  K  is  the  eddy  viscosity.  Therefore,  near  the 
equator,  bE  should  be  much  larger  than  the  Monin- 
Obukhov  length  scale,  which  is  nearly  20  m  in  the 
equatorial  Pacific  (Schneider  and  Muller  1990).  The 
upwelling  velocity  at  the  base  of  the  equatorial  Ekman 
layer  is  nearly  1  m  day-1  ( Wyrtki  1981 )  and  decreases 
upward.  Therefore,  the  vertical  velocity  in  the  OML 
under  weak  surface  wind  forcing  is  only  a  fraction  of 
the  vertical  velocity  at  the  base  of  the  Ekman  layer. 
The  value  of  0.1  m  day-1  might  be  a  reasonable  esti¬ 
mation.  Furthermore,  the  temperature  jump  at  the 
OML  base  Ts  -  T-h  is  on  the  order  of  3°C  or  less,  and 
the  net  ocean  surface  heating  at  the  equatorial  Pacific 
Qo  is  nearly  100  W  m-2  (Wyrtki  1981 ).  The  heating 
rate  caused  by  the  ocean  surface  warming  is  estimated 
by  Qo/ PwoCpwhw  ~  0.052°C  day-1.  The  cooling  rate 
caused  by  the  vertical  advection,  however,  is  estimated 
by  -w„h(Ts  -  T-h)/hw  «  -0.015°C  day-1,  which  is 
much  smaller  than  the  surface  heating  rate.  Therefore, 
it  is  reasonable  to  consider  the  OML  under  weak  sur¬ 
face  wind  forcing  as  a  warming  layer.  An  increase  in 
the  OML  thickness  hw  leads  to  a  decrease  of  dTJdt, 
which  implies  a  negative  correlation  between  dTs/dt 
and  hw. 

We  should  point  out  here  that  the  OML  model  only 
depicts  local  thermal  processes;  no  large-scale  ocean 
dynamical  processes  are  involved.  This  is  reasonable 
only  when  the  ocean  dynamics  is  assumed  to  come 
quickly  into  balance  on  an  interannual  time  scale  over 


which  the  SST  is  evolving  (Neelin  1991 ).  If  we  want 
to  know  both  the  slow  SST  mode  and  the  relatively 
fast  modified  ocean  wave  modes  or  if  we  want  to  im¬ 
prove  the  ENSO  prediction  model,  we  need  an  embed¬ 
ded  mixed  layer-ocean  circulation  model.  In  this  paper, 
(5),  (6),  and  (7)  are  base  equations  of  OML  Model. 

b.  WISHE  Model 

WISHE  is  a  multilayer  model,  in  which  the  tropo¬ 
sphere  is  divided  into  a  free  troposphere  and  a  subcloud 
layer.  Two  good  reasons  for  us  to  choose  WISHE  as 
the  atmospheric  part  are  1 )  realistic  thermodynamics 
(especially  the  treatment  of  cumulus  clouds  with  shal¬ 
low  updraft-downdraft  couplet)  for  the  investigation 
of  the  SST  mode  and  2 )  it  is  still  simple  for  numerical 
and  analytical  studies.  The  coordinate  system  is  chosen 
on  the  equatorial  /3  plane  with  x  in  zonal,  y  in  latitu¬ 
dinal,  and  z  in  vertical.  The  subgrid-scale  motions  are 
divided  into  the  deep  cumulus  convective  area  and  the 
quiescent  surrounding  environment.  In  addition  to  the 
deep  convective  updraft,  a  shallow  updraft-downdraft 
couplet  with  equal  updraft  and  downdraft  mass  fluxes 
is  also  included  (Fig.  1 ). 

The  subcloud  layer  of  thickness  h  obtains  moist  en¬ 
tropy  from  the  ocean  by  water  vapor  evaporation,  while 
it  loses  moist  entropy  to  the  lower  troposphere  both 
by  cloudy  downdraft  and  by  net  subsidence  within  the 
surrounding  environment.  The  free  troposphere  with 
thickness  H  obtains  moist  entropy  from  the  subcloud 


Stratosphere 


—  o’  —  - -  l-<r  - -  _  . 

Troposphere 


Fig.  1.  Schematic  presentation  of  the  AOSHE  model. 
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layer  by  shallow  updrafts  and  from  aloft  by  downward 
entropy  advection,  and  loses  by  net  radiative  cooling 
(Fig.  2).  Neglecting  the  local  rate  of  temperature 
change  inside  clouds,  the  following  set  of  equations  are 
given  by  (YE91 ): 

h  -j  (\ndeb)  =  C„|  V*|(ln0ev  -  \ndeh) 
at 

tryVVy(ln0e*  ln0„„ ) 

+  ( 1  -  o)wd(\n6eb  -  In 8em),  (8a) 

which  is  the  moist  entropy  balance  within  the  subcloud 
layer  with  thickness  h, 

d 

H—(\n8em)  =  -HR  +  <rsws(ln6eb  -  lndem) 
dt 

-  ( 1  -  a)wd(\n8eb  -  In Bem),  (8b) 

which  is  the  entropy  balance  within  the  free  troposphere 
with  thickness  H, 

g  J(ln0)  =  -iV2(l  -  <r)wd-gR,  (8c) 

which  is  the  conservation  of  dry  entropy  in  the  free 
troposphere, 

V-V* +  #-'[(  1  -  <T)wd+  awc]  =  0,  (8d) 

which  is  the  continuity  equation, 
d\h 

—jf  =  -V50*  -  k  X  pyYb  -  CDh-'  |V*| V*,  (8e) 


which  is  the  momentum  equation  on  an  equatorial  /3 
plane, 

5  In Beb  =  5  ln0,  (8f) 

*  m 

which  is  the  first  supplementary  relation  obtained  from 
assuming  the  equality  of  subcloud  layer  and  midtro- 
pospheric  moist  entropy  fluctuations  (YE91 ),  and 

Hb  =  ~cpaTbe5  In Beb  +  50,  (8g) 

which  is  the  second  supplementary  relation  obtained 
by  vertically  integrating  the  hydrostatic  relation. 

Here,  In Beb  is  the  moist  entropy  in  the  subcloud  layer; 
\ndem  the  pressure-weighted  mean  of  moist  entropy 
through  the  depth  of  the  troposphere;  ln0£,s.  the  saturated 
moist  entropy  of  the  ocean  surface;  8  the  potential 
temperature  in  the  subcloud  layer;  Yb  the  horizontal 
velocity  in  the  subcloud  layer;  and  50*  the  fluctuating 
component  of  geopotential  in  the  subcloud  layer.  The 
vertical  velocity  of  the  deep  cumulus  convective  area 
and  the  quiescent  surrounding  environment  are  wc  and 
wd,  respectively;  a  is  the  fractional  areal  coverage  of 
cumulus  convection.  In  addition  to  the  deep  convective 
updrafts,  a  shallow  updraft-downdraft  couplet  with 
equal  updraft-downdraft  mass  fluxes  is  considered.  The 
mass  flux  in  the  updraft  and  downdraft  are  denoted 
by  and  —  H  is  the  thickness  of  free  tropo¬ 
sphere;  Hm  is  taken  as  the  height  of  the  level  where  the 
equivalent  potential  temperature  0<,is  a  minimum  (see 
Fig.  2);  h  is  the  thickness  of  subcloud  layer;  N  the 
Brunt-Vaisala  frequency  of  the  troposphere;  R  the  net 
radiative  cooling  in  the  troposphere;  and  T*  and  Tm 


Strolosphere 


Fig.  2.  Entropy  exchange  processes,  sources,  and  sinks;  mean  entropy  distribution 
(after  Yano  and  Emanuel  1991 ). 
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are  the  moist  and  dry  adiabatic  lapse  rates;  cpa  is  the 
atmospheric  specific  heat  at  constant  pressure;  and  T 
and  Tb  are  temperatures  of  the  troposphere  and  the 
top  of  the  subcloud  layer,  respectively.  The  overbar 
(  _  )  in  (8g)  represents  a  mass- weighted  tropospheric 
average,  and 


€ 


Tb-T 

Tb 


(9) 


The  precipitation  efficiency  of  the  WISHE  model  is 
simply  defined  as  the  ratio  of  the  deep  updraft  mass 
flux  to  the  total  updraft  mass  flux: 


_  <7  Wc 

P  (awc  +  <rsws)  ' 
Solving  for  asws  gives 


( 1  -  tp) 

Ws  =  -  OWc 


(10) 


(11) 


3.  Ensemble  means  and  perturbations 

a.  Basic  state 

All  model  variables  are  decomposed  into  two  parts: 
Z  =  <Z>  +  Z', 

where  Z  represents  any  AOSHE  variable;  (Z  )  and  Z ' 
denote  ensemble  mean  and  perturbation,  respectively. 

A  basic  state  of  the  AOSHE  model  with  constant 
zonal  wind  U  is  evaluated  from  the  steady-state  solu¬ 
tions  for  Eqs.  (5),  (6),  and  (7)  of  the  OML  and  (8) 
of  the  WISHE  ( YE9 1 ) .  After  some  manipulations,  the 
mean-state  variables  are  given  by 


c2(Qo)  _  c,(uw*y 
pw0cpwST  ga(hw)8T  wK  e/’ 

(12a) 

(At)  =  v  +  ^  8T, 

Pw0^pw\'^w/  \fow/ 

(12b) 

-s: 

1 

£ 

II 

3? 

(12c) 

^  Hw*  ^  VCdP^o/ Pw0  1  ) 

( 1 2d) 

( 12e) 

(\ndem)=  <ln<U>-^^, 

& 

(12f) 

II 

( 12g) 

,  x  g(R) 

( 12h) 

<«>  =  c,/r'|(/|  in(<£>). 

(12i) 

Note  that  a  constant  zonal  wind  U  will  lead  to  a 
constant  zonal  gradient  of  OML  thickness  (Philander 
1990,  p.  108).  Therefore,  the  mean  OML  thickness 
(hf)  should  be  a  linear  function  of  x.  When  an  easterly 
basic  flow  appears  in  the  troposphere,  the  mean  OML 
thickness  (hw)  reduces  eastward.  In  the  present  study, 
however,  this  zonal  variation  is  absent.  It  is  acceptable 
when  we  concentrate  only  on  exploring  the  two-way 
correlation  between  temporal  change  in  SST  and 
mixed-layer  depth  perturbation  and  examining  its  ef¬ 
fects  on  the  behavior  of  the  unstable  mode  in  a  coupled 
air-sea  system.  One  should  be  very  careful,  however, 
when  applying  the  current  results  to  the  whole  Pacific 
basin. 

When  the  AOSHE  system  is  perturbed  from  its 
equilibrium  state,  the  air-ocean  feedback  mechanism 
between  the  WISHE  and  the  OML  either  causes  the 
perturbation  to  grow  (positive  feedback)  or  to  dampen 
(negative  feedback).  The  principal  purpose  here  is  to 
study  the  SST  mode  in  a  more  realistic  thermodynam¬ 
ical  coupled  system.  Hence,  the  energy  exchange  at  the 
air-ocean  interface  is  the  primary  focal  point.  There¬ 
fore,  we  shall  neglect  the  perturbation  of  the  temper¬ 
ature  advection  term  (more  related  to  wave  mode)  in 
OML  A't.  Hereafter,  we  omit  the  primes  for  pertur¬ 
bation  variables. 


b.  Criterion  for  different  surface  conditions 

After  the  decomposition  and  neglecting  of  the  high- 
order  terms,  the  continuity  of  momentum  flux  at  the 
ocean  surface  ( 3 )  becomes: 


Uw* 

^  uw*  ^ 


Ub_ 

\u\ 


sgn(  U), 


(13) 


which  indicates  that  in  the  mean  easterlies,  east  (or 
west)  wind  perturbation  leads  to  the  enhancement  (or 
reduction )  of  TKE  in  OML;  in  the  mean  westerlies, 
otherwise  is  expected.  The  criterion  for  the  different 
ocean  surface  conditions  becomes 


l+W\sgniU) 

—  (^w)^ 1  + 


<0o> 


1  + 


< hw ) 


>0  (14a) 


for  the  strong  surface  wind  forcing,  with  the  OML  in 
the  entrainment  regime;  and 


<k> 


—  1  + 


Qo 

(Qo) 


l  + 


<^> 


<0  (14b) 


for  the  weak  surface  wind  forcing,  with  the  OML  re¬ 
treating  to  the  Monin-Obukhov  length  scale  lw.  The 
surface  wind  and  heat  flux  perturbations  will  change 
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the  OML  regime  from  one  to  the  other.  At  the  ocean 
surface,  generally,  the  perturbation  heat  flux  Qo  is  much 
smaller  than  the  ensemble  mean  heat  flux  (Qo)-  The 
surface  wind  perturbation  ub  will  then  take  an  impor¬ 
tant  role  to  determine  the  OML  regime  since  the  back¬ 
ground  easterlies  are  rather  weak  ( |  C/j  ~  2  m  s_1, 
YE91 ). 

Four  important  processes  change  the  surface  heat 
flux:  SST-surface  heat  flux  feedback,  wind-surface  heat 
flux  feedback,  cloud-radiation  feedback,  and  evapo¬ 
ration-surface  heat  flux  feedback.  Generally,  SST- 
surface  heat  flux  feedback  and  wind-surface  heat  flux 
feedback  are  treated  as  the  most  important  thermal 
linkages  between  ocean  and  atmosphere  in  tropics.  The 
perturbation-surface  heat  flux  Q0  is  then  approximately 
estimated  by 

Qo  -PaoCpaCe[\U\  Ts  +  sgn  (U)ub((des)  -  (Beb))], 

(15a) 

where  Ce  is  the  ocean  surface  heat  exchange  coefficient; 
cpa  is  the  atmospheric  specific  heat  at  constant  pressure 
( YE9 1 ).  The  deep  water  temperature  perturbation  T_h 
is  determined  largely  by  the  physical  processes  in  the 
ocean  interior  and  is  assumed  negligible  in  the  present 
study.  According  to  (4)  the  entrainment  heat  flux  per¬ 
turbation  Q_h  is  then  computed  by 


dTs  1  _ 

<9/  —  ^ 

at  r  t 

[3(<2o>  +  PaoCpaCe\U\((6es)  -  (Beb))] 

Pw  o  cpw(hwy 

x~  sgn([/),  (17) 

where  rj'  is  thermal  dissipation  coefficient: 


rT  “ 


PwO^pw  <^> 
^PaOCpaQdl  0\ 


(18) 


2)  STRONG  SURFACE  WIND  FORCING 

During  strong  surface  wind  forcing,  the  OML  depth 
hw  satisfies  the  prognostic  equation  ( 7a) .  The  Eqs.  ( 5 ) , 
(6a),  and  (7a)  are  used  to  obtain  a  single  equation  for 
SST  perturbation.  If  the  vertical  velocity  at  the  OML 
base  w-h  is  assumed  to  be  determined  by  the  dynamics 
of  the  ocean  interior  and  not  to  be  directly  related  to 
the  surface  heat  exchange,  the  fluctuation  of  the  OML 
depth  for  a  strong  surface  wind  forcing  is  caused  by 
the  perturbation  of  entrainment  velocity  we, 


dhw 

dt 


(19) 


Q-h  =  PwoCpw((w-h)Ts  +  webT).  (15b) 

c.  Perturbations  in  the  OML  model 
1 )  Weak  surface  wind  forcing 


The  parameterization  of  we,  Eq.  ( 5 ) ,  indicates  that  the 
perturbation  entrainment  velocity  we  is  caused  by  the 
perturbations  of  the  SST  ( Ts),  the  OML  depth  hw,  the 
surface  water  friction  velocity  uw» ,  and  the  net  surface 
heat  flux  Q0, 


During  weak  surface  wind  forcing,  the  OML  depth 
is  taken  as  the  Monin-Obukhov  length  scale  (7b), 
which  is  the  algebraic  equation  for  the  three  variables: 
hw,  uw* ,  and  Q0.  The  fluctuation  of  the  OML  depth  is 
caused  by  the  perturbations  of  surface  heat  flux  and 
the  surface  wind  stress,  that  is, 

hw  _  3uw* _ Qo 

The  perturbation  heat  equation  becomes 


we 


C2Q0 

PwO^pw&T 


C\  ^  uw*  ^  hw 
ga(hw)dT  _(hw) 


-f^sgn  (tO 


(20) 


The  local  time  rate  change  of  perturbation  SST  is 
caused  by  perturbations  of  the  OML  thickness  hw,  of 
the  subcloud-layer  zonal  wind  (through  its  effect  on 
water  surface  friction  velocity),  and  of  surface  heat 
flux  Qo .  The  SST  perturbation  equation  for  strong  sur¬ 
face  wind  forcing  becomes 


dt 


(1  +C2) 
2 


tt'Ts  + 


(1  -C2)«2q)  +  2 (w-,)  g J  hw 
PwoCpw(hw'y  (hwy  _  (hwy 


C2(Qo )  +  (W-h)  dT 

pwo^ph'{^) 


+ 


(1  +  C2)paoCpaCe\U\((8es)  -  < 6eb » 
PwOCpw(  hy^y 


Ut 

\U\ 


sgn(C7).  (21) 


Substitution  of  ( 15)  and  (20)  into  (19)  leads  to  a  prognostic  equation  for  OML-depth  perturbation: 


dhw 

dt 


c2(Qo > 

PwoCpwbT 


hw 

<^> 


3|  C^Q^t+ {w-h) 

PwOCpwvl 


+ 


ClfiaO  CpqCfi  I  —  (Off,)) 

PwOCpw&P 


Ub 

\U\ 


sgn(f/).  (22) 
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A  relationship  between  SST  perturbation  Ts  and  surface 
wind  perturbation  C4  is  obtained  by  elimination  of 
OML-depth  perturbation  hw  from  (21 )  and  (22). 

4.  Air-ocean  feedback  mechanism 

The  tropical  subcloud  layer  is  usually  dominated  by 
mean  easterlies.  In  the  east  (west)  of  the  deep  convec¬ 
tive  region,  the  surface  wind  speed  is  enhanced  (re¬ 
duced);  therefore,  the  OML  depth  is  increased  (de¬ 
creased).  When  the  atmosphere  and  ocean  are  decou¬ 
pled  ( WISHE  model),  the  surface  evaporation  is  always 
enhanced  (reduced)  in  the  east  (west)  of  the  initial 
convective  region  by  the  WISHE  mechanism.  When 
the  atmosphere  and  ocean  are  coupled  (AOSHE 
model),  the  time  rate  change  of  surface  evaporation  is 
also  related  to  SST  (called  AOSHE  mechanism) 
through  the  first  term  in  the  right-hand  side  of  (8a). 
The  AOSHE  mechanism  is  different  from  the  WISHE 
mechanism  (e.g.,  enhanced  surface  evaporation  always 
in  the  east  of  the  convective  region)  in  such  a  way  that 
the  surface  evaporation  in  the  AOSHE  Model  also  de¬ 
pends  on  SST.  Generally,  the  increase  in  SST  augments 
the  surface  evaporation  rate.  The  effect  of  surface  wind 
and  buoyancy  forcing  on  SST  is  rather  complicated, 
however.  There  are  three  cases  available  for  the  tropical 
oceans. 

(i)  Weak  mean  easterlies  and  weak  initial  pertur¬ 
bations.  In  this  case,  the  total  surface  winds,  which  are 
the  summation  of  the  mean  easterlies  and  the  surface 
wind  perturbation  associated  with  the  initial  positive 
SST  fluctuation,  are  weak  such  that  OML  in  the  whole 
ocean  is  under  the  weak  surface  wind  forcing  ( Fig.  3a) . 


Therefore,  the  decrease  of  OML  depth  (hw  <  0)  leads 
to  the  augmentation  of  SST  perturbation  ( dTJdt  >  0) 
and  in  turn  to  a  production  of  the  atmospheric  con¬ 
vection  to  the  west  of  the  initial  convective  region  ( Fig. 
3a)  due  to  additional  moisture  convergence  and  surface 
evaporation  associated  with  the  augmentation  of  SST 
perturbation. 

(ii)  Strong  mean  easterlies  and  weak  initial  pertur¬ 
bation.  In  this  case,  the  total  surface  winds,  which  are 
the  summation  of  the  mean  easterlies  and  the  surface 
wind  perturbation  associated  with  the  initial  positive 
SST  fluctuation,  are  strong  such  that  OML  in  the  whole 
ocean  is  under  strong  surface  wind  forcing  (Fig.  3b). 
Therefore,  the  increase  of  OML  depth  ( hw>  0)  leads 
to  the  augmentation  of  SST  perturbation  ( 6Ts/dt  >  0) 
and  in  turn  to  a  production  of  the  atmospheric  con¬ 
vection  to  the  east  of  the  initial  convective  region  (Fig. 
3b)  due  to  additional  moisture  convergence  and  surface 
evaporation  associated  with  the  augmentation  of  SST 
perturbation. 

(iii)  Strong  initial  perturbation.  If  the  zonal  circu¬ 
lation  associated  with  the  initial  convection  is  strong 
enough  in  the  east  of  the  warm  water,  OML  is  under 
strong  surface  wind  forcing;  and  in  the  west  of  the  warm 
water,  OML  is  under  weak  surface  wind  forcing.  In 
the  east  of  the  warm  water  ( under  strong  surface  wind 
forcing),  the  increase  of  OML  depth  (hw>  0)  leads  to 
the  augmentation  of  SST  perturbation  ( 6Ts/dt  >  0) 
and  in  turn  to  a  production  of  the  atmospheric  con¬ 
vection  to  the  east  of  the  initial  convective  region  ( Fig. 
3c) .  In  the  west  of  the  warm  water  ( under  weak  surface 
wind  forcing),  the  decrease  of  OML  depth  (hw  <  0) 
leads  to  the  augmentation  of  SST  perturbation  ( dTs / 
dt  >  0)  and  in  turn  to  a  production  of  the  atmospheric 


Weak  Surface  Wind  Forcing  for  the  Whole  Area 


Fig.  3.  Atmospheric  convection  generated  by  the  air-ocean  feedback  mechanism  under  (a) 
weak  surface  wind  forcing,  (b)  strong  surface  wind  forcing,  and  (c)  strong  initial  perturbation. 
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convection  to  the  west  of  the  initial  convective  region 
(Fig.  3c).  Therefore,  when  the  initial  perturbation  is 
strong,  the  atmospheric  convection  will  be  generated 
in  both  east  and  west  sides  of  the  initial  convective 
area. 

The  unstable  modes  generated  by  this  air-ocean 
feedback  mechanism  are  called  AOSHE  modes,  which 
are  largely  determined  by  the  ocean  surface  condition. 
When  the  ocean  surface  is  under  weak  surface  wind 
forcing  this  mode  propagates  westward,  however,  when 
the  ocean  surface  is  under  strong  surface  wind  forcing 
this  mode  propagates  eastward.  Such  a  distinction  will 


largely  impact  upon  our  understanding  about  the  un¬ 
stable  equatorial  low-frequency  modes. 

5.  Nondimensional  system 

The  equations  are  nondimensionalized  following 
E87: 

x*  =  ax,  y*  =  a{/*A  l/4/?_l/2y,  z*  =  Hmz, 
h%  =  (hwyhw,  t *  =  a1/2A~l/1t, 
u*  =  al/2Al/2u,  v*  =  a~l/4A3/4l3~l/2v, 
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w*  =  HmAu2a~l/2w, 

d<t>*  =  icpaTbA(t>  =  aA<p,  8  ln0*  =  ATeb/t, 

8  In 0%  =  A Teb,  8  ln0?m  =  ATem, 

8  In**  «  -pr  =  ATa,  (23) 

\Is) 

where  the  asterisks  denote  the  dimensional  values.  Def¬ 
initions  and  typical  values  of  the  constants  are 

A  =  1cpaTbC)H~x  ln«0„>/<0rt» 

&  1.26  X  10-4  m  s~2, 

A  =  aH'lCe  \n((0es)/(0eb »  «  4.4  X  10“2, 

f  =  Td/Ym  «  1.7. 

Also  note  that 

A  =  (icpaTb/a)  A. 

Everything  in  (23)  is  the  same  as  in  YE91  except  that 
the  SST  perturbation  is  included. 

The  physical  values  of  the  model  constants  are  listed 
in  Tables  1  and  2.  From  these  values,  the  characteristic 
time  and  velocity  scales  are 

a1,2A~I/2  =  2.25  X  105  s  =  2.6  d, 

al/2Al/2  =  2$Ams~l.  (24) 

a.  Atmospheric  component 


Table  2.  Standard  values  for  oceanic  model  constants. 


PhO 

=  1.034  X  103  Kg  m-3 

Cpw 

=  4186  J  kg“'  K-' 

C, 

=  1.0 

c2 

=  0.2 

(hS 

=  50  m 

=  1  m  day-1 

6T 

=  3K 

(Ts) 

=  300  K 

m 

=  80  W  nr2 

ikub 

dvb 

dy 

F  ( 1  —  a)wd  +  awc  =  0, 

(25c) 

( D 

+  ocR)Teb  =  -fX(  1  -  a)wd. 

(25d) 

( 8D  + 

ae  +  aD)  Teb  -- 

=  sgn(  U)ub  +  aeTes 

+  X[ 

-o  - 

ep)awc 

“h  Cp(  1  &D Tem > 

(25e) 

( D  +  aD)Tem 

=  Ml 

-  ep)<rwc 

+  Wd 

-  (aR/t )]  Teb  -  cp\(  1  -  a)wd. 

(25f) 

where  the  terms  Tes,  <j>  are  used  for  the  coupling  with 
OML  and  the  stratosphere,  respectively.  The  effect  of 
the  stratosphere,  which  has  been  discussed  in  YE9 1 ,  is 
ignored  in  this  research,  that  is,  0  =  0  will  be  taken, 
corresponding  to  a  rigid  lid  at  the  tropopause.  Here, 
as  in  YE9 1 ,  the  model  solutions  are  assumed  to  have 
the  form 


From  the  foregoing,  a  nondimensional  set  of  equa¬ 
tions  for  the  atmospheric  part  is  given  by 

(£>  +  2F)ub  =  ik{Teb  -  0)  +  yvb,  (25a) 


(D  +  F)vb  =  P  \ 


d{  Teb  -  0) 
dy 


-  yub 


(25b) 


Table  1 .  Standard  values  for  atmospheric 
model  constants  (after  YE91). 


a  =  6.38  X  106  m 
0  =  2.3  X  10"8  s_1  km-1 

S  =  9.8m  s-2 
N2  =  10-4  s-2 
H  =  8  km 
Hm  =  5  km 

h  =  500  m 
CD  =  1  X  10~3 

C,  =  1.2  X  10-3 

_  ?  =  0.1 
In  (0J6*)  =  0.035 
\U\  =  2  m  s_l 
Cfa  =  1000  J  kg-‘  K-1 

Tmd  =  50  d 

Pa o  =  1-29  Kg  m-3 


radius  of  the  earth 
meridional  gradient  of  the 
Coriolis  parameter 
gravitational  acceleration 
Brunt-Vaisala  frequency 
thickness  of  the  troposphere 
level  of  minimum  Be  in  the 
troposphere 

thickness  of  the  subcloud  layer 
bulk  coefficient  of  the 
momentum  exchange  rate 
bulk  coefficient  of  the  entropy 
exchange 

thermodynamic  efficiency 
thermodynamic  disequilibrium 
magnitude  of  mean  zonal  wind 
air  specific  heat  at  constant 
pressure 

radiative  time  scale 
characteristic  air  density 


exp  (ikx  +  wt). 


and  also 

D^w  +  ikU,  P  =  pA~l/2ay2 

are  introduced.  Note  that  U  in  the  definition  of  D  is 
nondimensionalized  by  the  factor  Al/2ai/2.  The  pa¬ 
rameter  P  is  very  important  in  determining  the  dy¬ 
namical  feature  of  the  atmosphere.  If  P  is  large,  the 
zonal  wind  perturbation  is  in  geostrophic  balance;  if  P 
is  small,  the  meridional  wind  becomes  negligible,  and 
the  Kelvin-like  mode  dominates  the  atmosphere. 

The  nondimensional  parameters  that  appear  in 
(25a-f)  are  defined  with  approximate  magnitudes  by 
(see  YE9 1 ) : 

F=  CD\  U*\al/2A~l/2h~l  ~  0.9, 


8  -  h/H  ~  0.0625, 


X  -  N2Hm/g8  ~  1, 
aR  =  al/2A~l/2/TTad~  0.05, 
«,■  Ce\U*\ax/2A~'l2/H  ~  0.07, 


aD  ss 


al/2A~l/2Ce\U*\ 

H2N2 


0.03 


(26) 


These  parameters  are  kept  the  same  as  in  YE91  for 
easy  comparison  between  the  WISHE  and  AOSHE 
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models.  Detailed  discussion  on  the  physical  significance  The  roots  of  the  third-order  algebraic  equation 
of  these  parameters  can  be  found  in  YE9 1 . 

PAD)  =  0 


b.  Oceanic  component  with  a  weak  surface  wind 
forcing 

Nondimensionalization  of  the  perturbation  heat 
equation  under  weak  surface  wind  forcing  (17)  leads 
to 


(D  +  a<w)  -  ikU)Tes  =  -a<w)  sgn (U)ub,  (27) 
where 


„l/2  .-1/2 

- «  0.005  (28a) 

tt 


(„,)  _  [3<(2o>  +  PaoCpaCe\U*\((6es)  -  (6eb))]a 
"  Pw0cpw(hw)\U*\(Ts)A 

«  0.268  (28b) 


are  the  nondimensional  thermal  dissipation  coefficient 
and  the  air-ocean  coupling  coefficient,  respectively. 
The  superscript  (w)  indicates  the  weak  surface  wind 
forcing.  The  numerical  values  are  obtained  by  use  of 
the  model  constants  listed  in  Tables  1  and  2. 

Elimination  of  ub  and  Teb  from  the  atmospheric  mo¬ 
mentum  equations  (25a),  (25b),  and  the  SST  equation 
(27)  leads  to 


G(w)Tes 


■I  dVb  , 

ika3  — — I-  aiyvb 
dy 


.(w) 


sgn  (U), 


(29) 


which  shows  the  oceanic  heat  redistribution  caused  by 
weak  surface  winds.  Here  G(w),  at,  and  a3  are  fourth- 
order,  second-order,  and  first-order  polynomials  of  D, 

G{w)  =  (D  +  ajw)  -  ikU)P3(D) 

+  ika{uw)ae(D  +  aD )  sgn(C7), 

a\  =  (D  +  otp){8D  +  oid  +  cte) 

.  <xd<Xr  2  ,  D(D  +  aR ) 
+  -TT”  -  ao  + - S - , 


are  the  eigenvalues  of  the  WISHE  model  ( YE91 ). 

c.  Oceanic  component  with  a  strong  surface  wind 
forcing 

Nondimensionalization  of  a  set  of  OML  dynamical 
and  thermodynamical  equations  for  a  strong  surface 
wind  forcing  (21 )  and  (22)  leads  to 

[D  +  -  ikU]Tes  =  -ais)  sgn (U)ub  +  ahhw  (32) 


(D  +  p  -  ikU)h„ 


ST 


a  if  sgn  (U)ub  +  aTTes, 


(33) 


where 


als)  =  {l  +C2la'l2A  —  ~  0.003,  (34a) 


which  is  the  thermal  dissipation  coefficient,  and 
a 


a 


(■s)  = . 


|t/*l<rs>A 


CljQo )  +  ( W-h )  gj. 


+ 


(1  +  C2)PaoCpaCo\U\((des)  -  ( 6eb » 


0.18 


PwOCpw(  hKf 

(34b) 

which  is  the  air-ocean  coupling  coefficient.  The  su¬ 
perscript  (s)  indicates  the  strong  surface  wind  forcing. 
The  other  nondimensional  parameters  for  OML  in  (32) 
and  (33)  are  defined  by 

a'/2A~l/1 

i  - 

PwQCpW^hyf  (^hyf 


oth  = 


KTs) 

«  0.0059, 

olt  =  au2A~'/2 


C2A(Ts)tt 1  (w^h)A(Ts) 
28  T  ( hw)8T 


a3  =  D\{  \-tp),  (30) 

where  P3{D)  is  a  third-order  polynomial: 

P3(D)^BiD3  +  B2D2  +  B3D  +  BA,  (31a) 

where 

Bi  =  8  +  t~l, 

B2  =  B\ot£)  +  cte  +  otp)  +  2FB\  H - jj - , 

B3  =  aD(ae  +  aRt-')  +  2 F{B2  -  2 FBf 

+  ik  +  ( 1  —  ep)\k2 , 
B4  =  aD[2F(ae  +  orT-1)  +  ik].  (31b) 


-0.0343, 

m  -  a'/2A~l/{-^^hf  +  YrT-) 

\  PwOCpw \  “w/ O  •*  \*lw/  J 

«  0.018.  (34c) 

These  parameters  are  evaluated  by  a  set  of  physical 
values  for  the  model  constants  representing  the  general 
equatorial  Pacific  Ocean  conditions  shown  in  Table  2, 
where  (w_S£)  is  designated  the  vertical  velocity  at  the 
base  of  the  Ekman  layer.  From  Wyrtki’s  (1981)  esti¬ 
mation,  the  mean  upwelling  velocity  at  the  base  of  the 
Ekman  layer  (z  =  —8E)  is  nearly  1  m  day-1 .  Theoret¬ 
ically,  the  Ekman  depth  at  the  equator  is  infinity.  If  8E 
200  m  is  thought  to  be  a  reasonable  estimation  of 
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Ekman  layer  thickness,  the  mean  vertical  velocity  at 
the  OML  mean  depth  can  be  roughly  estimated  by 


/  \  w y 


0.25  m  day" 


Elimination  of  hw  from  (32)  and  (33)  leads  to 
(D  +  aiJ)  ~  ikU){D  +  p-  ikU)  Tes 

As)  TT\\  n  - L 


G{s)Tes  = 


.,  dvb 

ika-i  —  +  aiyvb 
dy 


a£'>Sgn(t/),  (36) 


which  shows  the  oceanic  heat  redistribution  by  the  sur¬ 
face  winds.  Here  G{s)  is  defined  by 

D  +  n  -  ^  ah  -  ikU^G(s) 

^(D  +  -  ikU)(D  +  n  -  ikU)Pi(D) 

+  ika(us)ae  sgn(  U)(D  +  aD) 

(Ts) 


X  I  D  ■+■  ft 


8T 


ah  -  ikU  .  (37) 


d.  Eigenvalue  problem 


where 


G(a2  -  ika3)  -  a^iD  +  2 F)au  sgn(U) 

= _ -  k2a3a4au  sgn(  U) 

(D  +  2F)a3(D  +  as-  ikU)P3(D) 

G[a2-P-'P3(D)(D  +  F)] 

2  (D  +  2 F)a3(D  +  as-  ikU)P3(D ) 

ai{D  +  as  -  ikU)P3(D) 


G[a2-P~'P3{D)(D  +  F)] 


strong  surface  wind  forcing;  a2 ,  a4  are  linear  functions 
of  D: 

a2  =  [Z>A(1  -  ep)ik  +  (D  +  aD)  sgn(t/)], 
a4  =  ae(D  +  ao). 


=  -ais)  sgn({7)^£>  +  u  -  <*h  ~  ikUjub 

+  ahaTTes.  (35) 

The  last  term  on  the  right-hand  side  of  (35),  ahaTTes, 
is  very  small  ( ahaT  ~  10-5)  and  will  be  neglected. 

Elimination  of  ub  and  Teb  from  the  atmospheric  mo¬ 
mentum  equations  (25a),  (25b),  and  the  SST  equation 
(35)  leads  to 


Let 


where 


/ \ 1/2 

V  =  *[231  y’  Vb  =  i)b^e 


(J»2/  2) 


3  = 


(1  +  4^2^3/^i)172  ' 


(40a) 


(40b) 


The  subcloud-layer  meridional  wind  equation  ( 38 )  be¬ 
comes 


d2vb 


+  3 


23  %?2 


-  jj2 \vb  =  0, 


dr\ 2  \ 

which  has  analytical  solutions 

vb  =  e-('2/2)Hn(v),  n  =  0,  1,  2, 
for 

3  -  23  §  =  2n  +  1, 


(41) 


(42) 


(43) 


Elimination  of  Tes  from  the  atmospheric  equation 
(25)  and  the  SST  equation  (29)  for  the  weak  surface 
forcing  or  SST  equation  (36)  for  strong  surface  wind 
forcing  becomes  a  single  equation  for  vb\ 


which  is  the  dispersion  relation.  Here  II n  (v)  is  the  Her- 
mite  polynomial  with  integer  order  n .  The  meridional 
wind  in  the  subcloud  layer  becomes 

vb  =  v„G i)e^2  =  e"'-3^2Hn(n) 

By  the  use  of  the  dispersion  relation  (43)  we  get 


vb  =  exp 


(nJTi  +  ;r2) 

2(2  n  +  1 


Uni  i 


23 


1/2 


y  ■  (44) 


1  —  <%3y2)vb  =  0,  (38) 


Thus,  in  order  to  get  a  well-behaved  solution, 
Re(— Wo,  for  n  =  0 


Re  «  +  I  >  for  n  =  2, 


(39) 


where  ( G ,  au,  as)  is  taken  as  ( G(w) ,  a(uw),  a(sw))  for 
weak  surface  wind  forcing  and  as  (G(s\  ais\  a(ss))  for 


(45) 


The  validity  of  equality  for  n  =  0  is  due  to  the  bound¬ 
ness  of  H0. 

The  primary  purpose  now  is  to  discuss  new  insight 
on  low-frequency  modes  brought  on  by  the  existence 
of  the  two-way  correlation  between  time  rate  of  change 
of  SST  and  OML  depth,  rather  than  to  do  a  mathe¬ 
matical  exercise.  Therefore,  only  the  zero-order  mode 
{n  =  0)  is  considered.  Substitution  of  (40b)  into  (43) 
for  n  =  0  leads  to 


(46) 
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which  is 

G2P\(D){D  +  F)2  -  PPi(D)(D  +  F) 

X  G[ika 3  +  aia4(Z)  +  2 F)au  sgn(£7) 

+  k2a2a^au  sgn (U)  +  a2G] 

+  P2  {Ga2[ika3  +  a1a4(D  +  2 F)au  sgn(t/) 

+  k2a2a^au  sgn(t/)]  -  (£>  +  2 F)ata3 

X(D  + as- ikU)Pl(D)}  =  0.  (47) 

For  a  small  value  of  P,  the  dispersion  relation  (47)  is 
simplified  as 

G2P\(D)(D  +  F)2  <=»  0,  (48) 

which  can  be  separated  into  three  algebraic  equations: 

P3(D)  =  0,  G(D)  =  0,  (D  +  F)  =  0, 

where  P3(D)  =  0  is  the  dispersion  relation  for  the 
WISHE  model;  D  +  F  =  0,  represents  the  damping 
modes  with  decay  rate  —  F  (not  of  interest  here);  and 

G(D)  =  0,  (49) 

which  is  the  dispersion  relation  for  new  modes  ap¬ 
pearing  only  in  the  AOSHE  model.  By  the  use  of  (39), 
this  solution  (49)  leads  to 

SP2  =  0, 

which  indicates  that  the  eigenvalue  solutions  are  well- 
behaved  solutions  [see  Eq.  (45)].  The  AOSHE  solution 
for  large  values  of  P  and  higher  orders  ( n  =  1,2,  •  •  • ) 
will  be  discussed  in  future  studies. 

6.  AOSHE  mode  under  weak  surface  wind  forcing 
(AOSHE  W-Mode) 

For  a  weak  surface  wind  forcing,  G  =  G(w),  and  the 
dispersion  relation  (49)  becomes 

(D  +  as-  ikU)P2(D) 

+  iky(w)  sgn(  U)(D  +  aD)  =  0,  (50) 

where 

7(M')  =  aLnV~0.02,  (51) 

which  comes  from  the  air-ocean  coupling.  The  value 
of  7 (M,)  is  estimated  by  (26)  and  (28b).  If  there  is  no 
air-sea  feedback,  y <lv)  =  0,  the  dispersion  relation  ( 50 ) 
is  reduced  to 

Ps(D)  =  0. 

Generally  speaking,  the  fourth-order  algebraic 
equation  (50)  has  four  roots.  For  a  given  set  of  non- 
dimensional  parameters,  only  two  of  the  four  roots  of 
(50)  have  positive  real  parts  representing  unstable 
modes.  The  remaining  two  or  more  roots  all  have  neg¬ 
ative  real  parts  representing  damped  modes,  which  are 
not  of  interest  here  and  are  therefore  neglected.  Among 


the  two  unstable  modes,  one  is  high  frequency  ( 30-60 
day),  and  the  other  is  low  frequency  (interannual). 
The  high-frequency  mode  is  the  modified  WISHE 
mode.  The  low-frequency  mode  is  called  the  AOSHE 
W-Mode,  that  is,  the  AOSHE  mode  under  weak  surface 
wind  forcing.  The  nondimensional  parameters  in  the 
dispersion  relation  (50),  y (lv),  a(sw),  X,  F,  and  ep,  will 
be  varied  in  the  present  study  in  order  to  investigate 
the  sensitivity  of  the  model.  The  physical  significance 
of  all  the  parameters  except  y (vv)  and  are  discussed 
in  YE91.  Here  yiw)  and  a*’1'*  denote  the  strength  of 
the  air-ocean  coupling  and  OML  thermal  dissipation. 

a.  Sensitivity  to  air-ocean  coupling  parameter  y (w) 

The  dispersion  relation  (50)  is  solved  for  the  eigen¬ 
value  D  for  a*M’)  =  0.005,  F  =  0.1,  ep  =  0.8,  X  =  1,  and 
for  five  different  values  of  7<M,):  0.01,  0.02,  0.03,  0.04, 
0.05.  For  7 {w)  =  0,  referring  to  the  air-ocean  decoupled 
case  (AOSHE  reduced  to  WISHE),  there  is  no  AOSHE 
W-Mode.  For  y <M,)  #  0,  the  AOSHE  W-Mode  is  gen¬ 
erated.  Figure  4  demonstrates  the  sensitive  dependence 
of  y(w)  of  the  AOSHE  W-Mode.  For  a  given  wave- 
number  k,  the  growth  rate  increases  monotonically 
with  7(lv). 

The  AOSHE  W-Mode  has  the  following  features. 

(i)  Largest-scale  mode  selection.  The  growth  rate 
Dr  has  its  maximum  at  the  lowest  wavenumber,  indi¬ 
cating  that  the  largest  scale  grows  most  rapidly.  The 
growth  rate  decreases  monotonically  as  k  increases.  For 
a  given  y(w),  there  exists  a  critical  wavenumber  kc. 
The  larger  scale  (k<kc)  AOSHE  W-Mode  is  a  growing 
mode;  however,  the  smaller  scale  (k>  kc)  AOSHE  W- 
Mode  is  a  damping  mode.  This  criterion  kc  increases 
monotonically  with  y (w):  kc  s=»  5  for  y iw)  =  0.01  to  kc 
sw  1 5  for  7 (K,)  =  0.05.  As  k  -*•  oo ,  the  growth  rate  tends 
to  a  negative  value  for  all  y (w),  which  shows  that  the 
smallest-scale  disturbance  damps  fastest  (a  distinct 
feature  from  the  WISHE  mode). 

( ii )  Slow  growth.  The  maximum  growth  rate  of  the 
AOSHE  W-Mode  (on  the  order  of  2.4  yr~'  for  y (w) 
=  0.02)  is  nearly  1.6%  of  that  of  the  WISHE  mode. 
The  growth  rate  increases  with  the  augmentation  of 
7  <w)>  which  depends  on  ( Q0 ) ,  Ce ,  A ,  < hw) ,  <  Ts) ,  H, 
and  A.  The  weaker  the  surface  wind  forcing,  the  shal¬ 
lower  the  OML,  and  in  turn  the  larger  the  value  of 
7 (w);  and  hence,  the  faster  growing  of  the  AOSHE  W- 
Mode. 

( iii )  Low  frequency.  The  AOSHE  W-Mode  is  a  low- 
frequency  mode.  The  period,  Pe  =  27t/|A:c|,  of  the 
unstable  mode  is  nearly  2-10  years  (Fig.  4b).  Such  a 
frequency  is  on  the  same  order  of  the  ENSO  frequency. 

(iv)  Slow  westward  propagation.  The  AOSHE  W- 
Mode  slowly  propagates  westward.  The  phase  speed  in 
all  cases  is  smaller  than  0.3  m  s-1  (Fig.  4c).  Similar  to 
the  growth  rate,  the  phase  speed  also  has  its  maximum 
at  the  lowest  wavenumber.  The  phase  speed  tends  to 
zero  as  the  wavenumber  k  tends  to  infinity.  The  phase 
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Fig.  4.  The  dependence  of  (a)  growth  rate,  (b)  period,  and  (c) 
phase  velocity  of  the  AOSHE  W-Mode  in  the  mean  easterlies  ( U* 
=  -2  m  s'1,  on  longitudinal  wavenumber  k  for  various  values  of 
the  air-ocean  coupling  parameter  y  <"  ).  The  other  important  param¬ 
eters  are  <*<"’>  =  0.005,  F  =  0.1,  tp  =  0.8,  and  X  =  1. 


speed  increases  with  the  augmentation  of  y{w).  The 
phase  speed  for  unstable  AOSHE  W-Mode  ( 10-30 
cms"1)  shows  a  quite  good  agreement  with  the  ob¬ 
served  westward  phase  speed  (20  cm  s'1 )  of  the  equa¬ 
torial  Pacific  14°C  isotherm-depth  perturbation  (Mey¬ 
ers  1979)  and  with  the  observed  westward  phase  speed 
(30-50  cm  s_1 )  of  the  surface  wind  and  ocean  signals 
(Michum  and  Lukas  1987). 


The  solutions  show  that  there  is  almost  no  effect  of 
air-ocean  coupling  on  the  high-frequency  mode  ( 30- 
60  day),  which  means  that  only  the  WISHE  mecha¬ 
nism  is  responsible  for  the  generation  of  the  intrasea- 
sonal  variation. 

b.  Sensitivity  to  a{w> 

The  dispersion  relation  (50)  is  solved  for  the  eigen¬ 
value  D  for  7(m,)  =  0.02,  F  =  0.1,  ep  =  0.8,  X  =  1,  and 
for  five  different  values  of  a(s"’>:  0,  0.002,  0.005,  0.01, 
0.02.  Figure  5  indicates  the  sensitivity  of  the  AOSHE 
W-Mode  growth  rate  on  the  thermal  dissipation  coef¬ 
ficient  aiM). 

(i)  For  a  given  wavenumber  k,  the  growth  rate  de¬ 
creases  monotonically  with  For  alw>  =  0  (no 
thermal  dissipation  in  OML),  the  growth  rate  is  always 
positive  and  tends  to  zero  as  k  >— »  oo . 

( ii )  When  a  i w)  >  0,  a  wavenumber  criterion  kc  exists 
such  that  the  larger  scale  (k  <  kc)  AOSHE  W-Mode  is 
growing,  and  the  smaller  scale  ( k  >  kc)  AOSHE  W- 
Mode  is  damping.  This  criterion  kc  decreases  mono¬ 
tonically  with  the  increasing  ajM,)  from  kc  15  for 
a{sw)  =  0.002,  to  kc  ss  2  for  a\w)  =  0.02. 

Furthermore,  a  very  slight  effect  of  ai*1  on  the 
AOSHE  W-Mode  phase  speed  was  found.  In  the  largest 
scale  k  <  3,  the  augmentation  of  the  thermal  dissipation 
a[w)  increases  the  phase  speed  by  a  few  percent.  For  k 
>  3,  there  is  no  effect  of  a(sw)  on  the  AOSHE  W-Mode 
phase  speed. 

c.  Sensitivity  to  the  mean  easterlies  V* 

The  annual  mean  easterlies  in  the  central  and  eastern 
Pacific  are  about  5  m  s'1,  while  in  the  western  Pacific 
they  are  about  2msH.In  order  to  address  the  model 
sensitivity  with  respect  to  U*,  we  solve  the  dispersion 


Fig.  5.  The  dependence  of  the  growth  rate  of  the  AOSHE  W- 
Mode  in  the  mean  easterlies  ( U*  =  —  2  m  s'1 )  on  longitudinal  wave- 
number  k  for  various  values  of  the  OML  thermal  dissipation  a\w>. 
The  other  important  parameters- are  yw  =  0.02,  F  =  0.1,  cp  =  0.8, 
and  X  =  1.0. 
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relation  (50)  for  F  =  0.1,  y(w)  =  0.02,  a<w)  =  0.005, 
tp  =  0.8,  X  =  1 ,  and  for  five  different  values  of  U*:  —1 
m  s_1 ,  -2  m  s“‘ ,  -3  m  s-1 ,  -4  m  s-1 ,  and  -5ms'1. 

Figure  6  shows  the  dependence  of  growth  rate  and 
phase  speed  of  the  AOSHE  W-Mode  on  k  for  different 
values  of  U*.  It  is  noticed  that  the  growth  rate  (Fig. 
6a)  for  the  whole  wavenumber  domain  and  the  phase 
speed  (Fig.  6b)  for  k  >  3  are  insensitive  to  the  mean 
wind;  however,  for  k  <  3,  the  westward  phase  speed 
increases  with  the  increasing  U*.  Not  much  variation 
in  both  growth  rate  and  phase  speed  for  different  values 
of  U*  indicates  that  the  AOSHE  mechanism  is  valid 
for  the  whole  ocean  basin  (western,  central,  and  east¬ 
ern)  when  OML  is  in  the  shallowing  regime. 

7.  AOSHE  mode  under  strong  surface  wind  forcing 
(AOSHE  S-Mode) 

For  a  strong  surface  wind  forcing,  G  =  G(s),  and  the 
dispersion  relation  (49)  becomes 

(D  +  «<s)  -  ikU)(D  +  n-ikU)P3(D) 

=  -ikyU)  sgn(U)(D  +  aD)(D  +  g  -  ye  -  ikU), 

(52) 


Fig  .  6.  The  dependence  of  the  (a)  growth  rate  and  (b)  phase  velocity 
of  the  AOSHE  W-Mode  on  longitudinal  wavenumber  k  for  various 
values  of  the  mean  easterlies  U*.  The  other  important  parameters 
are  y (B,)  =  0.02,  a<“')  =  0.005,  F  =  0.1,  X  =  1,  and  tp  =  0.8. 


where 

7<s>  =  a\pae  —  0.01,  7e-a*^-~0.1,  (53) 

which  are  two  important  parameters  representing  the 
ocean-atmosphere  surface  coupling  and  the  OML  en¬ 
trainment  process.  Here,  the  values  of  y (i)  and  ye  are 
estimated  by  the  use  of  (26),  (34b),  (34c),  and  Table 
2.  It  is  noticed  that  if  there  is  no  air-sea  feedback,  y (i) 
=  0,  the  dispersion  relation  (52)  is  reduced  to  P3{D) 
=  0,  that  is,  the  AOSHE  model  reduces  to  the  WISHE 
model.  The  physical  significance  of  y  e  denotes  the 
strength  of  the  entrainment  mixing.  The  larger  the 
value  of  y  e,  the  stronger  the  entrainment  process  in 
the  OML. 

Generally  speaking,  the  fifth-order  algebraic  equa¬ 
tion  (52)  has  five  roots.  For  a  given  set  of  nondimen- 
sional  parameters,  at  most  two  of  the  five  roots  of  ( 52 ) 
have  positive  real  parts  representing  unstable  modes. 
Among  them,  one  is  similar  to  the  WISHE  solution 
(i.e.,  high-frequency  mode),  and  the  other  is  a  new 
solution  only  appearing  in  the  AOSHE  model  and  rep¬ 
resenting  a  low-frequency  unstable  mode  called 
AOSHE  S-Mode.  The  other  three  roots  all  have  neg¬ 
ative  real  parts  representing  dampening  modes,  which 
are  not  of  interest  here  and  are  therefore  neglected. 
The  five  nondimensional  parameters  in  the  dispersion 
relation  (52),  7 <s),  ye,  a*s),  F,  tp,  will  be  varied  in  the 
present  study  in  order  to  investigate  the  sensitivity  of 
the  model. 

a.  Sensitivity  to  air-ocean  coupling  parameter  y (s> 

The  dispersion  relation  (52)  is  solved  for  the  eigen¬ 
value  D  for  ye  =  0.1,  =  0.006,  F  =  0.1,  tp  =  0.8, 

X  =  1,  and  for  five  different  values  of7(s):  0.005,  0.01, 
0.02,  0.03,  and  0.04. 

The  growth  rate  and  phase  speed  remain  unchanged 
for  all  7 (s)  for  the  high-frequency  (30-60  day)  mode. 
This  indicates  that  the  WISHE  mechanism  is  respon¬ 
sible  for  the  generation  of  the  high-frequency  ( 30-60 
day)  mode.  The  air-sea  coupling  is  not  important  for 
such  a  mode. 

Figure  7  demonstrates  that  the  sensitive  dependence 
of  7 <J)  actually  leads  to  the  generation  of  the  low-fre¬ 
quency  (interannual)  AOSHE  S-Mode.  For  y =  0, 
referring  to  the  air-ocean  decoupled  case,  the  AOSHE 
mode  does  not  appear;  for  y {s)  +  0,  the  AOSHE  S- 
Mode  does  appear. 

The  AOSHE  S-Mode  has  the  following  features  (Fig. 
7): 

(i)  Length- scale  mode  selection.  The  growth  rate  Dr 
has  its  maximum  in  the  wavenumber  domain.  The 
maximum  growth  rate  shifts  toward  the  lower  wave- 
number  side  as  y <s)  increases,  that  is,  the  wavenumber 
with  maximum  growth  rate,  /cmax ,  changes  from  kmax 
5  for  7 <s)  =  0.005,  to  /cmax  ^  4  for  y (s)  =  0.03. 
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Fig.  7.  Dependence  of  the  low-frequency  AOSHE  S-Mode  on  the 
coupling  parameter  7 (J>:  (a)  growth  rate,  (b)  period,  and  (c)  phase 
speed.  The  other  important  parameters  are  ye  =  0. 1,  ai'1  =  0.006, 
F  =  0.1,  e„  =  0.8,  and  X  =  1. 


( ii )  Slow  growth.  The  maximum  growth  rate  of  the 
AOSHE  S-Mode  is  0.46  yr~’  for  y (I)  =  0.01,  which  is. 
much  smaller  than  the  AOSHE  W-Mode. 

(iii)  A -type  effect  on  growth  rate.  It  is  interesting  to  . 
see  the  A-type  effect  of  the  air-ocean  coupling  coeffi¬ 
cient  y (s)  on  the  growth  rate  of  the  AOSHE  S-Mode: 
the  growth  rate  generally  increases  with  increasing  y (s) 
for  0.005  <  y(s)  <  0.01,  and  then  decreases  with  in¬ 
creasing  y <s>  for  0.7  <  y(s)  <  1.5. 


(iv)  Low  frequency.  The  AOSHE  S-Mode  is  a  low- 
frequency  mode.  The  period,  Pe  =  2ir/\kc\,  of  the 
most  unstable  mode  is  nearly  1-3  years  (Fig.  7b).  The 
period  increases  with  wavenumber  k,  and  decreases 
with  the  air-ocean  coupling  parameter  y(j).  Such  a 
frequency  is  on  the  same  order  of  the  ENSO  frequency. 

(v)  Slow  eastward  propagation.  The  AOSHE  S- 
Mode  slowly  propagates  eastward  (Fig.  7c).  The  phase 
speed  is  in  the  range  of  0- 1 .25  m  s ~ 1 .  For  a  given  y (i), 
the  phase  speed  decreases  monotonically  with  k,  and 
tends  to  zero  as  k  -*■  oo.  For  y is)  =  0.01,  the  phase 
speed  decreases  from  85  cm  s-1  for  k  =  1  to  12  cm  s-1 
for  k  =  5.  For  a  given  k,  the  phase  speed  increases 
monotonically  with  y(s). 

b.  Sensitivity  to  entrainment  parameter  y  e 

The  dispersion  relation  ( 52 )  is  solved  for  the  eigen¬ 
value  D  for  7<J)  =  0.01,  a]s)  =  0.006,  F  =  0.1,  ep 
=  0.8,  X  =  1,  and  for  five  different  values  of  ye:  0.01, 
0.05,0.1,0.25,0.5. 

Figure  8  demonstrates  that  the  sensitive  dependence 
of  ye  actually  contributes  to  the  generation  of  the  low- 
frequency  (interannual)  AOSHE  S-Mode. 

(i)  For  small  values  of  ye  (ye  «£  0.05 )  representing 
a  weak  entrainment,  the  growth  rate  for  the  AOSHE 
S-Mode  is  always  negative.  This  means  that  the  en- 


Fig.  8.  Dependence  of  the  low-frequency  AOSHE  S-Mode  on  the 
entrainment  parameter  ye:  (a)  growth  rate  and  (b)  phase  speed.  The 
other  important  parameters  are  y <J)  =  0.01,  =  0.006,  F  =  0.1, 

tp  =  0.8,  and  X  =  1. 
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trainment  mixing  in  the  OML  should  be  sufficiently 
strong  to  generate  unstable  modes.  For  any  wavenum¬ 
ber  k,  the  growth  rate  increases  monotonically  with  y  e 
(Fig.  8a).  The  maximum  growth  rate  of  the  AOSHE 
S-Mode  increases  0.46  yr_1  for  y  e  —  0.01  to  2.5  yr-1 
for  ye  =  0.5. 

( ii )  The  AOSHE  S-Mode  propagates  eastward  ( Fig. 
8b).  For  a  given  ye,  the  phase  speed  decreases  mono¬ 
tonically  with  k,  and  tends  to  zero  as  k  -*■  oo .  For  a 
given  k,  the  phase  speed  increases  monotonically  with 
ye.  As  y  e  =  0.01,  the  phase  speed  is  nearly  zero;  how¬ 
ever,  as  y  e  =  0.1,  the  phase  speed  reduces  from  85 
cm  s“‘  at  k  =  1  to  5.8  cm  s-1  at  k  -  8. 

c.  Sensitivity  to  thermal  dissipation  coefficient  a[s> 

The  dispersion  relation  ( 52 )  is  solved  for  the  eigen¬ 
value  D  for  y  w  =  0.01,  ye  =  0.1,  F  =  0.1,  tp  =  0.8,  X 
=  1,  and  for  five  different  values  of  a  *s):  0, 0.003, 0.006, 
0.01,  0.02.  Figure  9  indicates  the  sensitivity  of  the 
AOSHE  S-Mode  growth  rate  on  the  thermal  dissipation 
coefficient  a(ss). 

(i)  For  a  given  wavenumber  k,  the  growth  rate  de¬ 
creases  monotonically  with  a{ss).  For  a|.s)  =  0  (no  ther¬ 
mal  dissipation  in  OML),  the  growth  rate  is  the  largest. 
For  =  0.02  (approximately  3.3  times  the  regular 
thermal  dissipation  coefficient),  there  is  no  unstable 
mode  at  all  for  the  whole  wavenumber  domain. 

(ii)  The  unstable  AOSHE  S-Mode  has  length-scale 
mode  selection.  The  maximum  growth  rate  is  located  at 
k  5  for  all  a  The  size  of  the  wavenumber  interval 
for  the  unstable  modes  increases  as  a{ss)  decreases. 

A  very  slight  effect  of  a  [s)  on  the  AOSHE  S-Mode 
phase  speed  was  found.  In  the  largest  scale  k  <  3,  the 
augmentation  of  the  thermal  dissipation  a(sw)  increases 
the  phase  speed  by  a  few  percent.  In  the  smaller  scale 
(k  >  3),  there  is  no  effect  of  ctlss)  on  the  AOSHE  S- 
Mode  phase  speed. 


Fig.  9.  Dependence  of  the  low-frequericy  AOSHE  S-Mode  growth 
rate  on  the  thermal  dissipation  parameter  a1/1 .  The  other  important 
parameters  are  7W  =  0.01,  ye  =  0.1,  F  =  0.1,  e,  =  0.8,  and  X  =  1. 


Fig.  10.  The  dependence  of  the  growth  rate  (a),  and  phase  velocity 
(b)  of  the  AOSHE  S-Mode  on  longitudinal  wavenumber  k  for  various 
values  of  the  mean  easterlies  U*.  The  other  important  parameters  are 
7(J>  =  0.01,  7„  =  0.1,  a<s)  =  0.006,  t,  =  0.8,  X  =  1,  and  F  =  0.1. 


d.  Sensitivity  to  mean  easterlies  U* 

The  dispersion  relation  (52)  is  solved  for  y  (5)  =  0.01, 
y e  =  0.1,  =  0.006,  F  —  0.1,  X  =  1,  ep  =  0.8,  and 

for  five  different  values  of  U*\  -lms"',-2m  s~' , 
-3  m  s-1,  —4  m  s~‘,  and  —5  m  s-1. 

Figure  10  shows  the  dependence  of  growth  rate  and 
phase  speed  of  the  AOSHE  S-Mode  on  k  for  different 
values  of  U*.  It  is  noticed  that  the  growth  rate  (Fig. 
10a)  for  a  given  wavenumber  decreases  with  increasing 
values  of  the  mean  easterlies  U*  and  that  this  reduction 
is  enhanced  as  wavenumber  k  increases.  It  is  also  no¬ 
ticed  that  the  phase  speed  remains  almost  the  same  for 
different  values  of  the  mean  easterlies  U*.  The  reduc¬ 
tion  of  growth  rate  with  the  increasing  U*  implies  that 
the  AOSHE  mechanism  weakens  from  the  western  Pa¬ 
cific  to  the  central  and  eastern  Pacific  when  OML  is 
in  the  entrainment  regime. 

8.  Different  model  sensitivities  between  AOSHE  W 
and  S  modes 

The  different  types  of  sensitivity  of  AOSHE  W  and 
S  modes  on  the  model  parameters  (sections  6  and  7) 
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are  outlined  in  Tables  3  and  4.  Here,  the  upward  arrow 
f  ( or  downward  j )  indicates  the  increase  ( or  decrease ) 
of  the  growth  rate  (Table  3)  or  phase  speed  (Table  4) 
with  the  increase  of  parameters.  The  parameters  (y, 
a5)  are  ( y(w) ,  a]K’))  for  the  AOSHE  W-Mode,  and 
(7 (w>,  a<">)  for  the  AOSHE  S-Mode. 

Three  important  atmospheric  parameters  A,  F,  and 
tp  have  totally  different  effects  on  the  growth  rate  and 
phase  speed  of  the  two  AOSHE  modes.  The  increase 
of  F  generally  represents  the  enhancement  of  drag  at 
the  ocean  surface,  CD\U*\  [see  Eq.  (26)],  which 
strengthens  the  TKE  generation  in  OML.  More  TKE, 
generated  by  the  surface  wind  stress,  in  turn  increases 
the  OML  depth,  which  leads  to  dTs/dt  <  0  in  the 
AOSHE  W-Mode,  and  dTJdt  >  0  in  the  AOSHE  S- 
Mode.  The  two  distinct  effects  of  precipitation  effi¬ 
ciency  tp  on  the  AOSHE  model,  augmentation  for 
AOSHE  W-Mode  and  attenuation  for  AOSHE  S- 
Mode,  are  also  due  to  the  two  different  kinds  of  OML- 
thickness  SST  feedback.  The  precipitation  brings 
buoyancy  to  OML  and  tends  to  reduce  the  OML  depth, 
regardless  of  whether  OML  is  under  weak  or  strong 
surface  wind  forcing.  In  other  words,  the  precipitation 
will  induce  a  negative  OML  depth  fluctuation  ( hw 
<  0),  which  causes  the  two  different  tendencies  for  the 
SST  perturbations:  dTs/dt  >  0  for  a  weak  surface  wind 
forcing,  and  dTs/dt  <  0  for  a  strong  surface  wind  forc¬ 
ing.  These  two  distinct  types  of  SST  response  have  two 
different  effects  on  the  atmospheric  convection:  en¬ 
hancement  when  dTs/dt>  0  (weak  surface  wind  forc¬ 
ing)  and  attenuation  when  dTs/dt  <  0  (strong  wind 
forcing  surface). 

For  the  two  AOSHE  modes,  the  air-ocean  coupling 
coefficient  y  has  different  effects  on  growth  rate  and 
the  same  effect  on  phase  speed;  however,  the  OML 
thermal  dissipation  coefficient  as  has  the  same  effect 
on  growth  rate  and  different  effects  on  phase  speed. 

9.  Summary 

(i)  There  are  two  distinct  air-ocean  feedback 
mechanisms  available  for  the  two  different  tropical 
ocean-surface  conditions.  As  the  ocean  surface  is  under 
weak  wind  forcing,  the  time-rate  change  of  SST  is  neg¬ 
atively  correlated  with  OML  depth  hw;  however,  as  the 
ocean  surface  is  under  strong  wind  forcing,  the  time- 
rate  change  of  SST  is  positively  correlated  with  OML 
depth  hw.  Such  a  difference  leads  to  the  generation  of 
two  different  low-frequency  (interannual)  modes  in  the 


Table  3.  Different  effects  of  model  parameters  on  the  growth  rate 
of  the  two  AOSHE  modes. 
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Table  4.  Different  effects  of  model  parameters  on  the  phase  speed 
of  the  two  AOSHE  modes. 
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mean  easterlies;  the  eastward-propagating  AOSHE  S- 
mode  (under  strong  surface  wind  forcing)  and  the 
westward-propagating  AOSHE  W-mode  ( under  weak 
surface  wind  forcing).  This  air-ocean  feedback  mech¬ 
anism  is  called  an  AOSHE  mechanism,  which  might 
play  an  important  role  in  the  ENSO  events. 

( ii )  The  AOSHE  model  predicts  the  generation  of 
a  westward-  (or  eastward)  propagating  low-frequency 
(interannual)  mode  under  weak  (or  strong)  surface 
wind  forcing.  By  the  use  of  typical  values  of  the  pa¬ 
rameters  for  the  tropical  Pacific,  the  westward-propa¬ 
gating  AOSHE  W-mode  shows  the  following  features: 
largest-scale  mode  selection,  faster  growth  (maximum 
growth  rate  is  on  the  order  of  (2.2  yr_1 ),  low  frequency 
(2-10  year  period  for  the  unstable  mode),  slow  west¬ 
ward  propagation  (phase  speed  on  the  order  of  -20 
cm  s'1);  however,  the  eastward-propagating  AOSHE 
S-Mode  indicates  the  following  properties:  wavenum¬ 
ber  3-5  mode  selection,  slower  growth  (maximum 
growth  rate  on  the  order  of  0.5  yr-1),  low  frequency 
( 1-3  year  period  for  the  most  unstable  mode),  faster 
eastward  propagation  (phase  speed  10-100  cm  s_l). 

(iii)  The  AOSHE  mechanism  is  valid  for  the  whole 
ocean  basin  (western,  central,  and  eastern)  when  OML 
is  in  the  shallowing  regime;  however,  the  AOSHE 
mechanism  weakens  from  the  western  Pacific  to  the 
central  and  eastern  Pacific  when  OML  is  in  the  en¬ 
trainment  regime. 

(iv)  Besides  the  AOSHE  W-  and  S-modes,  the 
AOSHE  model  preserves  the  modified  WISHE  modes. 
The  air-sea  feedback  effect  on  the  modified  WISHE 
mode  is  negligible.  This  indicates  that  the  AOSHE 
mechanism  and  the  WISHE  mechanism  are  apparently 
not  interfering  with  each  other.  The  WISHE  mecha¬ 
nism  is  responsible  for  the  generation  of  the  high-fre¬ 
quency  mode  (30-60  day)  in  the  atmosphere  with 
smallest-scale  mode  selection;  and  the  AOSHE  mech¬ 
anism  is  responsible  for  the  production  of  the  low-fre¬ 
quency  mode  (interannual)  with  either  largest-scale 
mode  selection  for  the  AOSHE  W-Mode,  or  wave- 
number  3-5  selection  for  the  AOSHE  S-Mode. 

( v)  In  the  present  study,  the  advection  perturbations 
are  neglected.  The  system  depicts  local  thermal  pro¬ 
cesses  only.  There  are  no  large-scale  ocean  dynamic 
processes  involved.  The  ocean  waves  are  filtered  out 
in  the  current  AOSHE  model.  Furthermore,  OML  with 
salinity  effect  is  not  included  here.  Combination  of  the 
ENSO  wave  theory  and  the  AOSHE  mechanism  with 
consideration  of  salinity  effect  will  provide  a  complete 
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picture  of  dynamical  and  thermodynamical  interaction 
between  equatorial  troposphere  and  the  ocean. 

(vi)  The  AOSHE  model  can  also  be  served  as  a 
linkage  between  the  ocean  and  the  stratosphere  since 
the  stratosphere  was  already  included  into  the  WISHE 
Model  in  YE9 1 .  The  inclusion  of  the  stratosphere  into 
the  AOSHE  model  might  lead  to  a  possible  linkage  of 
the  low-frequency  oscillations  in  the  ocean  and  in 
the  stratosphere,  such  as  the  quasi-biennial  oscilla¬ 
tion  (QBO). 

(vii)  The  advantage  of  the  equality  assumption  of 
subcloud  layer  and  midtropospheric  moist  entropy 
fluctuations  (8f),  is  to  simplify  the  model  without  loss 
of  major  physical  processes;  however,  this  assumption 
becomes  too  strong  if  one  wants  to  develop  a  numerical 
ENSO  prediction  model  including  the  AOSHE  mech¬ 
anism. 
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